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ABSTRACT 

Since their discovery, cosmic crystalline silicates have presented several chal- 
lenges to understanding dust formation and evolution. The mid-infrared spec- 
trum of IRAS 17495—2534, a highly obscured oxygen-rich asymptotic giant 
branch (AGB) star, is the only source observed to date which exhibits a clear 
crystalline silicate absorption feature. This provides an unprecedented opportu- 
nity to test competing hypotheses for dust formation. Observed spectral features 
suggest that both amorphous and crystalline dust is dominated by forsterite 
(Mg2Si04) rather than enstatite (MgSiOs) or other silicate compositions. We 
confirm that high mass-loss rates should produce more crystalline material, and 
show why this should be dominated by forsterite. The presence of Mg2Si04 glass 
suggests that another factor (possibly C/0) is critical in determining astromin- 
eralogy. Correlation between crystallinity, mass-loss rate and initial stellar mass 
suggests that only the most massive AGB stars contribute significant quantities 
of crystalline material to the interstellar medium, resolving the conundrum of its 
low crystallinity. 

Subject headings: stars: AGB and post-AGB — stars: individual (IRAS 17495- 
2534) — (stars:) circumstellar matter — (ISM:) dust, extinction — infrared: 
stars 



1. Introduction 

One of the most exciting recent developments in astronomy w as the discovery of crys 



talline silicate Stardust by the Infrared Space Observatory (ISO; Ide Graauw et al.l Il996l ) 
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Crystalline silicates were i nitially discovered a round evolved intermediate mass stars at far- 
infrared (IR) wavelengths ([Waters et al. 119961) . but have since been detected in young stellar 



objects (IWaelkens et a. 



Galaxies (ISpoon et al. 



19961) ■ comets (ICrovisier et al.lll996l ). and Ultra Luminous Infrared 
2006^ In the inter stellar medium ( I SM) silicate grains i n crystalline 



form are rare (~1% crystalline by mass; iMin et al.l 120071 : iKemper et al.l 12004 ) . However, 
Asymptotic Giant Branch (AGB) stars and their successors (pre-planet ary and planetary 



nebulae) are the major contributors of material to the ISM (iKwok |2004J) . It has been sug- 
gested that AGB stars could contain ■^10% crystalline silicates (IKemper et al.l 120041 ). so 



it is vital to understand the formation of dust around AGB stars to resolve this apparent 
discrepancy. 

AGB stars are lumin ous cool giants, with h igh mass- loss rates (10^^ < M < IQ-'^Mq/jt), 
which increase over time (Ivan Loon et al.ll2005l ). They are highly evolved descendants of low- 
to intermediate-mass (0.8-8 M0) stars. Mass loss produces a circumstellar envelope, where 
dust forms. The order in which different species condense from the outflowing gas depends 
on the physical conditions within the envelope. Once M is very high the circumstellar shell 
becomes optically thick, even at IR wavelengths, and the silicate emission features absorb 
themselves. Th e dust mineralogy is d i ctated by the dust conde nsation sequence and depends 
strongly on M (jPijkstra et al.ll2005l ). iBlommaert et al.l (120071 ) showed that the mineralogy 
depends only on the age of the star for a given initial stellar mass. 

To date, crystalline silicate features have been seen clearly only in emission. It has 
been proposed that crystalline silicate absorption at llfim contributes to the broadening 
of the classic lO/zm amorphous silicate absorption feature associated with OH/IR stars 
JSvlvester et al.lll999l : IVanhoUebekel I2OO7I ). However, IRAS 17495-2534 (hereafter 117495) 
is the first object to show a distinct crystalline 11/zm absorption feature. 

117 495 is located in t he Galactic Plane (gal. coord. = 003.6844+00.3880) at a distance of 
~4kpc (ILoup et al.lll993l ). about half way to the Galactic Center. The IRAS LRS spectrum 
clearly exhibits an ll.l^m absorption feature superposed on the classic lOyum amorphous 
silicate feature (Fig. [1]). This is the first clear crystalline silicate absorption feature seen to 
date in any object. 

M for this object is estimated to b e ~ 2 x 10"^ < M < 5 x lO'^'Mo/yr based on 
mid-IR ([25]- [12]) color (ILoup et al.lll993l ) and the range of optical depths (15 < rio^m < 45) 
consistent with o ur model results (see §[3]). Both M and expansion velocity (vexp, 16km/s; 



Loup et al.lll993l ) are at the high end of normal for AGB stars. Most extreme 0-rich AGB 
stars are OH/IR stars, exhibiting OH-maser emission. 117495 is one of a handful of visibly 
obscu red 0-rich AGB stars not exhibiting OH-maser emission (dubbed color-mimics; [Lewis 
1994[ ). However, the spectra of other color mimics more closely resemble those of OH/IR 
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stars, i.e., when observed, the putative 11/im crystaUine sihcate feature is merely a hump 
superposed on the regular amorphous silicate absorption feature. This is also true for OH/IR 
stars with similar mass-loss rates in the Galact ic Bulge (GB). Even the GB source with the 
highest M (3 X lO-^Mg/yr; IRAS 17276-2846; IVMihoUebekd |20o3) shows the 11/xm feature 
as a hump, rather than the distinct 11/im feature seen in Fig. [H 

While forste rite (crystalline Mg2Si04) is expected to have a peak near 1 1. Sum (e.g. 
Koike et al.ll2003h . the peak can shift to shorter wavelengths. In the laboratory, iJager et al.l 



( 1l998l ) found that synthetic forsterite peaks closer to 11.2/ i ni, and matche s the p osition of the 
hump observed in GB OH/IR stars ( VanhoUebekd 120071 ) . iFabian et al.l ( 200ll ) investigated 
grain shape effects. Their mass absorption coefficients for ellipsoidal forsterite grains are 
included in Fig . [ Tl and demonstrate that the crystalline forsterite feature can peak at ll.l//m. 



Boersma et al. 



(EoOS ) att ributed an ll.lum em ission feature in the spectrum of a Herbig 



Ae/Be star to forsterite. iTamanai et al.l (120061 ) found that, for non-embedded free-flying 
particles of forsterite, the feature actually peaks at 11.06/im. Composition, grain-shape and 
grain-size can signiflcantly shift the position of this feature. 



2. Crystallinity & Mass-loss rates 



Crystalline silicates were flrst observed around evolved stars w ith very high M, leading 
to the inference that crystal formation requires such conditions (e.g. ICami et al.lll998l ). How- 
ever, the diagnostic infrared features of crystalline grains can rernain h idden by the strong 
amorphous silicate features until M 10~^MQ/yr (IKemper et al.l 1200 ll ). 



Dust condensation temp erature (T^ust) are depressed at low gas pressures (ILodders fc Fegley 



19991 : iGail fc Sedlmayrlll999l ). Fig. [2] shows the calculated pressure -temperature (P-T) con- 
ditions in the dust-forming zone around 0-rich AGB stars (see 



19991) 



Diikstra et al. 


2005. for 


n models ( 


Lodders & Fedev 


(Richet et al. 


1993 


) and 



MgSiOs (IWilding et al.l 120041 ). Grains forming below Tg will be amorphous, while grain s 



forming much above this temperature should crystallize rapidly (ISogawa fc Kozasal Il999l ). 
For M 2 X 10"^ Mg/yr, Tdust for both Mg2Si04 and MgSiOg fall below their respective 
Tg, thus precluding crystal formation. Above this threshold M, dust grains should form 
at temperatures where crystalline structures would form either directly, or by annealing of 
initially amorphous grains. Once dust forms, the radial temperature proflle steepens, be- 
cause of absorption of starlight by dust grains. Including this factor, we can use the outflow 
velocity to calculate the timescale for dust to cool from T^ust to Tg. For M > lO~^M0/yr, 
the timescale is ~6-12 months, providing ample time for annealing and crystallization. This 



-4- 



suggests that crystalline silicates should be more common in high M stars and the observed 
M-c rystallinity correlation i s not merely due to masking of the crystalline spectral features 
(c.f. ISogawa fc Kozasal Il999l ). In addition, for stars with M^10~^ Mq /jt forsterite forms 



at a higher temperature than enstatite (crystalline MgSiOa) 
these two species increases with M (see §|1]). 



and the difference in TH„<=t for 



3. Radiative transfer (RT) modeling 



We used the 1-D radiative transfer (RT) program DUSTY (jivezic fc Elitzuiill995l ). to 
determine the gross miner alogy associated with t he mid-IR absorption features. RT modeling 
is inherently degenerate (jivezic fc Elitzurl 119971 ) and hence models are not unique. Conse- 
quently, we use the modeling only to constrain rio/im and to determine which minerals are 
inconsistent with the observed spectral features. A detailed analysis of the model results 
will be presented elsewhere. 

At pressures relevant to 117495, Tdust Ril200K for Mg2Si04 and MgSiOs. Models with 
Tdust < 800 K could not match the observations because the near-IR becomes too bright, and 
in the mid-IR the spectrum slopes the wrong way. One important result of the modeling is 
that pyroxenes can be excluded from the mineralogy in 117495. The underlying amorphous 
silicate feature is unu sually narrow arid pea ks closer to 10/xm rather than 9.7/im, consistent 
with Mg2Si04 glass (lOssenkopf et al.lll992l . and references therein ) . Models incorporating 
amorphous silicate of pyroxene cornpositi on (e.g. iJager et al.lll994l ) or other more silica-rich 
compositions (e.g. lOssenkopf et al.l Il992l ) are inconsistent with the observed spectrum, as 
the peak position of the feature gets too blue, and the FWHM gets too wide. Increasi ng the 
iron content in olivine has a similar effect on the peak position (lOssenkopf et al.lll992l ). Our 
models suggest that the crystalline component is also dominated by forsterite with enstatite 
and other crystalline compositions limited to less than a few percent of the dust mass. 
Therefore, both amorphous and crystalline dust is dominated by Mg2Si04. This observed 
forsteritic mineralogy needs to be understood in terms of dust formation mechanisms. 



4. Competing dust formation mechanisms 



There are effectively three competing dust formatio n mechanisms for circurn stellar envi- 
ronments: (i) thermodynamic equilibrium condensation (ILodders &: Feglevlll999f): fill fo rma- 



tion of chaotic solids in a supersaturated gas followed by annealing (IStencel et al.lll990[): (Hi) 



formation of seed nuclei in a supersaturated gas, followed by mantle growth (iGail &: Sedlmayr 
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19991 ). The latter should follow thermodjTiamic equilibrium as long as density is high enough 
for gas-grain reactions to occur. 



Several observational stu dies support the thermodynamic condensation sequence (iDijkstra et al. 



20051 : iBlommaert et al.ll2007l ). which is consistent with both (i) and (iii). In mechanism (ii), 
chaotic grains fo rm with the bulk co mposition of the gas, and then anneal if the temperature 
is high enough (IStencel et al.lll990l ). This mechanism predicts that at low C/0 ratios, the 
dust grains would comprise a mixture of olivine, pyroxene and silica, rather than be domi- 
nated by olivine alone. At high C/0 ratios, Al-0 bonds are predicted to form preferentially, 
leading to dust dor ainated by oxides rather than silicates. Thes e predictions are inconsistent 
with observations (iDijkstra et al.ll2005l : IBlommaert et al.ll2007l ). 



In the seed-mantle model, Mg2Si04 is expected to use more refractory condensates as 
nucleation centers (e.g. AI2O3; see Fig. Whether the transformation into MgSiOs by 
gas-grain reaction occurs is then dependent on gas density. Condensation and annealing 
experiments show that the resulting dust i s dominated by forsterite; enstatite is present only 
in trace amounts (e.g. iDay fc Donnlll978l : iNuth fc Donrull982l : iHallenbeck et al.lll998l ). In 
laboratory experiments, the density is always high compared to circumstellar shells, so they 
may only be relevant to the highest M. 

The modeled mineralogy is consistent with an abrupt decrease in density between the 
formation temperatures for forsterite and enstatite. Once ~20% of the Si atoms are incor- 
porated into dust, the opacity is such that radiation pressure accelerates the circumstellar 
mater ial, resulting in a preci pitous density decrease, and inhibiting further gas-grain reac- 
tions (IGail fc Sedlmayiill999l ). At high M, forsterite is stable at a higher temperature than 
enstatite. Formation of forsterite grains provides the impetus for radiation-driven accelera- 
tion and thus reduces the efficacy of enstatite formation, leaving gas that is depleted in Mg 
and thus more silica-rich. Without a surface on which to nucleate, silicates will tend to form 



amorphous solids rather than crystals (iTsukamoto et al.ll200ll ). thus the seed nuclei are of 
the utmost importance to forming crystalline grains. Moreover, the dispersion due to radia- 
tion pressure may lead to seedless amorphous dust formation in the low density gas, which 
should be slightly more silica-rich than MgSiOa due to the depletion of Mg into forsterite 
grains. The glass composit ion should tend t o wards Mg2Si308, which is one of the metastable 
compositions predicted by iRietmeijer et al.l (119991 ). 



It is possible that the appearance of the clear crystalline absorption feature in 117495 
is related not to increased crystallinity, but rather to less maski ng of the feature as a result 
of similar crystalline and amorphous silicate compositions (see iKemper et al.ll200ll ). This 
still leaves the question of why this source is different - how can this source have both its 
crystalline and amorphous dust components dominated by Mg2Si04, given that the bulk 
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composition of the gas is closer to MgSiOs? 



5. Factors influencing crystallinity and mineralogy 



The physical factors which could give rise to the extremely high crystal content include: 
(1) M; (2) fexp; (3) metalhcity; (4) C/0 ratio. M for 117495 is high, but similar to OH/IR 
stars and other color mimics, suggesting that M alone cannot explain the high degree of 
crystallinity exhibited by this source. A low v^^^ would increase the pressure at a given 
temperature and thus promote grain formation/growth at higher temperatures. However, 
for 117495 is relatively high. Increasing metallicity increases the partial pressure of 



the dust forming elements, and thus increases the dust formation temperature at a given 
total gas pressure. The metallicity of the source is unknown, although its location in the 
Galactic plane suggests that it may have higher than solar metallicity. If metallicity were 
the determining cause, we would expect to see similar spectra in the GB OH/IR stars. The 
fact that even the high M GB stars do not exhibit such a distinct crystalline silicate feature 
suggests that a combination of M and metallicity would also fail to increase the crystal 
fraction. 

Having eliminated M, fexp and metallicity as likely causes of higher crystallinity, we are 
left with C/0 ratio as the most likely culprit. The condensation scheme suggested in § H] 
is relevant to solar C/0 (~0.48). Oxygen-rich AGB stars have C/0 ratios between 0.4 and 
0.8, where the exact value depends on both the initial mass and age of the star; increasing 
this ratio decreases the availability of oxygen atoms with which to make dust. 

The Mg/Si ratio in the circumstellar outflow from which the dust forms is close to unity. 
The silica-poor composition of the amorphous dust requir es sequestration of sora e silicon 
into a non-silicate phase to increase the Mg/Si ratio (c.f. iFerrarotti &: Gaill l200ll ). While 
the formation of gehlenit e fCa^Al^SiO?) before the Mg-rich silicates is predicted in m any 
condensation models (e.g 
references therein), Si/Ca 



Lodders fc Feglevlll999l : iGail fc Sedlmavi]ll999l : lTielenslll99d . and 



^16, and thus gehlenite is not an effective sink for silicon atoms. 



Other options would be keeping the Si in gas phase SiO; or condensation of Si into a 
metallic (or other featureless) phase (e.g. FeSi). The predicted path by which MgSiOs forms 
is a reaction of Mg2Si04 with SiO gas, requiring an additional oxygen atom. Consequently, 
high C/0 may limit the conversion of Mg2Si04 into MgSiOs and pres erve a silica-poor 
mine r alogy. The conderisation of FeSi is predicted for C/0 ratios >0.65 (IFerrarotti fc Gail 
2OO2I : iLodders fc FeglevI I1999I). As C/0 approaches unity, Tdust for the Mg-silicates de- 
creases rapidly (ILodders fc Fegleyl 119971 ) and Si is more apt to be incorporated into FeSi 
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( iFerrarotti fc Gaill |2002| ). For C/Q > 0.85, the Mg-silicate dust formation temperature 
drops to below that of FeSi (IFerrarotti &: Gai]||2002l ). which results in sequestration of Si and 
increases the Mg/Si ratio of the gas from which the sihcates form. Therefore, there are at 
least two mechanisms by which higher C/0 ratios may be able to explain the composition 
of 117495. 

Based on the preceding arguments we suggest that a high C/0 ratio is the cause of 
the unusual mineralogy of 117495. Our hypothesis predicts several observable correlations: 
(1) between the strength of SiO absorption and C/0 ratio in low M stars; (2) between 
SiO absorption and the relative abundances of crystalline forsterite and enstatite in high M 
stars; and (3) between C/0 and the peak positions/FWHM of the classic 10/xm (emission) 
feature. As discussed in § E] the peak position and FWHM of the lO^m feature depend on 
the composition of the amorphous silicate. 



6. The effect of source morphology 



Moist er et al.l (119991 . |2002| ) found a correlation between AGB/post-AGB star crystallinity 



and source morphology , with "disk-l i ke" so urces having a higher degree of crystallinity 
than spherical sources. iMolster et al.l (119991 ) suggested that this correlation is due to low- 
temperature crystallization in a long-lived disk. However, the inference of a disk around 
highly crystalline sources comes from imaging and the general shape of the SED, which 
cannot distinguish between a Keplerian disk and a toroidal outflow. There is evidence that 
the last stages of mass loss from AGB stars become increasing toroidal, an d that this ef- 
fect is stronger the more massive the progenitor star (jDijkstra fc Speckll2006l . and references 
therein). In this case, the density /pressure in the dust-forming regions would be enhanced 
further by the concentration of material into the star's equatorial region, leading to con- 
ditions more apt to form crystalline structures. This in turn would correlate degree of 
crystallinity with initial mass of the star. Only massive AGB stars would produce significant 
amounts of crystalline material, and these are significantly rarer than low mass AGB stars. 
This resolves the conundrum of the ISM (having low crystallinity) being enriched by high 
crystallinity sources (see §[1]). 



7. Conclusions 



We have presented the first crystalline silicate absorption feature observed to date. The 
spectrum of 117495 is enigmatic and suggests that its dust is dominated by Mg2Si04 in 
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both the crystalline and amorphous phases. We have confirmed that high mass-loss rates 
should produce more crystalline material than low mass-loss rates; and that the crystalline 
mineralogy should be increasingly dominated by forsterite as mass-loss rates increase. The 
most likely factor controlling both crystallinity and mineralogy is the C/0 ratio. We suggest 
that the correlation between crystallinity, mass-loss rate and initial stellar mass mitigates 
the problem of the very different crystal fractions observed in some AGB stars and in the 
ISM. Only the rarer, higher mass AGB stars contribute a significant amount of crystalline 
material to the ISM. 

We are grateful to H.-P. Gail, whose comments significantly improved this paper. This 
work is supported by NSF CAREER AST-0642991 (for A.K.S.) and NSF CAREER EAR- 
0748411 (for A.G.W.). 
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Fig. 1. — The spectrum of IRAS 17495— 2534: y-axis is flux (AFa) in Wm~^. Main panel: 
log-log plot of the full spectral energy distribution (SED) including the IRAS LRS spec- 
trum, the near-IR JHK photometry from 2MASS, 3.6/im photometry form the Spitzer 
GLIMPSE survey, and far-IR photometry points from IRAS. Inset: linear-scale close-up 
on the continuum-divided IRAS LRS spectrum, with crystalline features indicated by ar- 
rows, y-axis is normalized flux. Grey lines are calculated mass absorption c oefficients for 



contin uous distributions of ellipsoidal (CDE) crystalline forsterite grains from iFabian et al. 



(I2OOII ). CDEl has equal probability of all shapes; CDE2 has quadratic weighting, such that 
near-spherical particles are most probable. Dotted line at 9.7/um is included to demonstrate 
the redness of the peak of the absorption feature. 
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Fig. 2. — P-T space for dust condensation around 0-rich AGB stars: Solid and dashed lines 
indicate Tjust for a given pressure from thermody namic equilibrium calculations (relevant 
compositions are labeled) (ILodders fc Fegleylll999l ). The stability temperature of AI2O3 is 
dot-dashed where extrapolated. For all M values, AI2O3 forms at a significantly higher 
temperature than the silicates, and thus can form a seed nucleus. Dotted lines indicate the 
P-T paths for the outflowing gas for a range of M (indicated in MQ/yr) as calculated from 
M and t'exp- Thick grey horizontal lines indicate Tg for Mg2Si04 and MgSiOs. Grey shaded 
area shows the P-T space relevant to 117495. 



